Mechanisms in the size segregation of a binary granular mixture 



O 

o 

(N 
o 

Q 



o 



Matthias SchroterIB Stephan Uhich, Jennifer Kreft, Jack B. Swift, and Harry L. Swinney 
Center for Nonlinear Dynamics and Department of Physics, 
The University of Texas at Austin, Austin, Texas 78712, USA 
(Dated: February 6, 2008) 

A granular mixture of particles of two sizes that is shaken vertically will in most cases segregate. 
If the larger particles accumulate at the top of the sample, this is called the Brazil-nut effect (BNE); 
if they accumulate at the bottom, the reverse Brazil-nut effect (RBNE). While this process is of 
great industrial importance in the handling of bulk solids, it is not well understood. In recent years 
ten different mechanisms have been suggested to explain when each type of segregation is observed. 
However, the dependence of the mechanisms on driving conditions and material parameters and 
hence their relative importance is largely unknown. In this paper we present experiments and 
simulations where both types of particles are made from the same material and shaken under low air 
pressure, which reduces the number of mechanisms to be considered to seven. We observe both BNE 
and RBNE by varying systematically the driving frequency and amplitude, diameter ratio, ratio of 
total volume of small to large particles, and overall sample volume. All our results can be explained 
by a combination of three mechanisms: a geometrical mechanism called void filling, transport of 
particles in sidewall-driven convection rolls, and thermal diff'usion, a mechanism predicted by kinetic 
theory. 

PACS numbers: 45.70.Mg, 05.20.Dd 
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I. INTRODUCTION 



A variety of mechanisms have been proposed to de- 
scribe the separation of particles of two sizes in a mixture 
of vertically shaken particles. Section II discusses pro- 
posed mechanisms, including void filling, static compres- 
sive force, convection, condensation, thermal diffusion, 
non-equipartition of energy, interstitial gas forcing, fric- 
tion, and buoyancy. Our aim is to determine which are 
the important segregation mechanisms for diverse condi- 
tions. 

We study a mixture of two sizes of brass spheres, using 
image analysis to count the number of large particles vis- 
ible at the top and bottom surfaces of the sample after 
shaking. These results are then compared with molecular 
dynamics simulations, which yield quantities unaccessi- 
ble in experiment. Comparing our results for a range 
of control parameters with various proposed segregation 
mechanisms reveals which mechanisms are important in 
our system. The control parameters studied are: fre- 
quency / and amplitude A of the sinusoidal shaking, the 
diameter ratio of the spheres dj^/ds (the subscripts L and 
S refer to large and small particles), the total layer depth 
htoted, and Vl/Vs, the ratio of the total volume of large to 
small particles in the container. The last parameter has 
not been examined in previous experiments, which were 
either performed with about equal numbers of particles 
of each size [l|, 0, H, 0, S 0], or in the case called the 
intruder limit, where there is one lar ge sphe re in a bed 
of smaller particlesM, i, i, [M [HI, [Illll 113, M, M, 113, 



This paper is organized as follows: different proposed 
segregation mechanisms are reviewed in Sec. |TT1 Section 
mil discusses the interplay of these mechanisms in two 
recent models [H, . The experimental setup and the 
details of our simulation are described in Sec. lIVI Section 
|V] discusses how the experimental and simulation results 
depend on the driving parameters, and Sec. I VII discusses 
the dependence on di,/ds, Vl/Vs and /itotai- We conclude 
in Sec. Em 



II. SEGREGATION MECHANISMS 



To reduce the number of segregation mechanisms to be 
considered we use large and small particles of the same 
material (brass) and therefore the same density and me- 
chanical properties. This exclude s segregati on effects due 
to Archimedean buoyancy H, 23 
different frictional properties 28 



251 1261, |27|,|28| or due to 
29|. As the interstitial 
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fluid was also found to have a strong effect on segregation 
d, [H, [13, nil, [23, 1131 , we reduced the pressure in our ex- 
periment to less than 8 Pa, well below the 130 Pa where 
the effect of air on the segregation vanishes '^I'l. This 
step also assures comparability with our molecular dy- 
namics simulations where the only interactions between 
particles are collisions. The remaining seven mechanisms 
are reviewed below, and their predicted dependence on 
dh/ds, Vl/Vs and fetotai is summarized in Table HI For a 



recent review see also 31 
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A. Geometrical mechanism 

1. Void filling 

Void filling is a local geometric effect that was first 
clearly identified in Monte Carlo simulations [s^. Dur- 
ing the expansion phase of the shaking cycle all particles 
move upwards, but during the compaction phase a large 
particle has a smaller probability of finding a suitable 
void in the layer below it than a small particle. This 
leads to the large particles accumulating at the top of 
the sample. 

As pointed out in (sij , void filling works only for weak 
excitation; for strong driving the voids become larger and 
the probability that they accommodate a large particle 
becomes similar to that of a small particle. The void 
filling effect becomes stronger with increasing diameter 
ratio di^/ds, as shown by the increase in the rise velocity 
of an intruder 13]. For d^/ds ~ 3 the intruder seems to 
undergo a transition from intermittent to continuous as- 
cent, which indicates a change in the microscopic details 
of the mechanism Q . Comparing in a thought exper- 
iment the situation of a single large intruder in a bed of 
small particles with a mixture which contains predomi- 
nantly large particle hints that increasing Vl/Vs will de- 
crease the strength of void filling. In contrast, the overall 
layer height /itotai should not influence this mechanism. 

There exists no theoretical model for void flUing, but 
Trujillo et al. [25] identify in their kinetic theory based 
analysis of the gravitational force a "dynamic tensile" 
force; for an intruder this upward force increases with 
[d-L/d^Y. This term vanishes in the limit of a highly 
expanded mixture. They compare the dynamic tensile 
force with the void filling mechanism. 



2. Static compressive force 

Trujillo et al. ^25^ find that in the presence of gravity, 
any size disparity will give rise to another term, which 
they name static compressive force. The static compres- 
sive force will always lead to the RBNE; in the intruder 
limit it is proportional to (dL/'^s)'^ — 1- 



B. Convection 

Friction at the container walls can induce convection 
rolls in vertically shaken mixtures: during the upward ac- 
celeration the mixture gets compacted and shear forces 
induced by the side walls extend through the whole sam- 
ple. During the subsequent downward motion the mix- 
ture is more expanded and consequentially those parti- 
cles adjacent to the walls experience stronger shear forces 
than those in the center of the container. Both phases 
combined give rise to a convection roll going downwards 
at the sidewalk and upwards in the center. 



The relevant geometrical control parameter seems to 
be the ratio of the size of the downstream layer to the di- 
ameter of the large particles. If the large particles are too 
large to be entrained by the downstream, they become 
trapped at the top surface and consequently convection 
leads to BNE [ifl, [O, El. However, if the size of the 
downward moving regions is big enough to accommodate 
larg e particles, these particles perform convective cycles 
[la Q ini, m, [11], which can lead to mixing Q. 

Magnetic resonance imaging indicates that the size of 
the downward moving layer is independent of total layer 
depth and shaking acceleration [ll[. However, increasing 
^totai might reduce the efficiency of both mechanisms as 
the strength of the convection roll vanishes at the bottom 
[3^. Increasing the ratio of larg e to small particles Vl/Vs 
in a gedanken-experiment [3l| increases the size of the 
downstream layer and therefore suppresses the BNE and 
fosters mixing. 



C. Kinetic energy driven mechanisms 

For driven granular media one can define a kinetic 
granular temperature T which corresponds to the aver- 
age kinetic energy of the random motion of the particles 



^=y((«-(^'))'), (1) 

where m is the particle mass and v their velocity. Nu- 
merical studies [H, 113, HI] and experiments [H, 53, El 
have shown that T in a vertically shaken granular sam- 
ple depends not only on the position in the sample and 
the phase of the shaking cycle, but also on the shaking 
acceleration and the number of particles in the container. 

In binary granular mixtures the equipartition of en- 
ergy breaks down [i^, [4ll, [Z^, [4^ ; in general the heavier 
particles have a higher granular temperature than the 
lighter ones. Non-equipartition increases with decreasing 
restitution coefficient and increasing mass difference 44] . 

Two different types of segregation mechanisms driven 
by the granular temperature of the sample have been 
suggested, based either on a phenonienological model of 
solid-liquid phase transitions [2^ or on analyzing mo- 
mentum balance equations derived from kinetic theory 

[Ulillii,!!!, 113,141,14^. 



1. Condensation 

A mechanism suggested by Hong et al. [23] is based 
on the observation that as the shaking acceleration is in- 
creased, a solid-like sample of monodisperse particles will 
become vibro-fluidized. Assuming that the granular tem- 
perature is spatially homogeneous throughout the sam- 
ple, they postulate the existence of a critical temperature 
T'^ above which all particles are fluidized. For T smaller 
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than T'^ an increasing fraction of the particles should con- 
dense at the bottom of the container. By equating the 
average kinetic energy of the particles with the potential 
energy of a particle in the top layer, they find T'^ to be 
proportional to pd^h, where p is the material density and 
h is the height of the particular layer. 

Hong et al. argued that in a binary sample the 
shaking parameters can be adjusted such that T is be- 
tween the critical temperatures of the two species. Then 
one kind stays fluidized while the other condenses at the 
bottom. Assuming cquipartition and that the values of 
T'^ are not influenced by the presence of the other species, 
they identify the ratio of the two T'^'s as the control pa- 
rameter econd, 
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where the height ratio hi^/hg of the particles filled in 
separately is equal to Vi,/Vs in a container of constant 
width. 

For econd < 1, which corresponds to < Tg, the 
small particles are still condensed and sink to the bot- 
tom while the large particles are completely fluidized; 
the system shows BNE. For econd > 1 the large particles 
will condense at the bottom and the RENE takes place. 
However, these effects will only occur in an intermedi- 
ate driving regime where T is in between the two critical 
temperatures; otherwise there will be no segregation as 
a result of condensation. 



2. Mechanisms predicted by kinetic theory 

Kinetic theory of binary granular mixtures is an ac- 
tive area of research with competing predictions. As ki- 
netic theory is limited to a regime where binary collisions 
prevail, these mechanisms are only applicable for dilute 
mixtures and hence large shaking amplitudes. 

Diffusion in a granular temperature gradient. It 
was first shown by Hsiau and Hunt that in a gradient 
of the granular temperature the heavier particles move to 
places with lower T and thereby expel the lighter parti- 
cles. Strong vertical temperature gradients in vertically 
shaken samples have been found in simulation [sgI . IstI . [38j 
and experiment [s^, |43| . Due to the complicated spatio- 
temporal development of the T field, it is not possible to 
decide a priori what type of segregation this effect will 
lead. 

Segregation due to non-equipartition. The 

derivation in (45| assumes equipartition. Recent sim- 
ulations by Galvin et al. I44l| have shown that non- 
equipartition leads to additional driving forces that are 
comparable in strength to the forces in the calculations 
assuming equipartition. 

An explicit solution for the thermal diffusion co- 
efficients in the case of a single intruder and non- 
equipartition was given by Brey et al. [isj . who found 



TABLE I: Effect of an increase of the diameter ratio, volume 
ratio and total layer depth on the segregation mechanism dis- 
cussed in Sec|lll The symbols 0, and - correspond to an 
increase, Constance or decrease of the effect in the second col- 
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the ratio of mean square velocities. 



(3) 



to be the control parameter. For particles with same 
density, they predict cj) < 1, which leads to RBNE; (j) 
decreases with increasing diameter ratio. 

Trujillo et al. [25j identify in their hydrodynamic model 
a so-called pseudo-thermal buoyancy force, which is pro- 
portional to the ratio of the granular temperatures of the 
two species. They calculate this temperature ratio using 
the kinetic theory results of Barrat and Trizac [13], where 
Tl is always larger than Tg so that the pseudo-thermal 
buoyancy leads to BNE. Barrat and Trizac also find that 
non-equipartition increases strongly with di^/ds and de- 
creases weakly with V^/Vs. 



III. COMBINATION OF MECHANISMS 

To predict the state of segregation in an experiment 
the different mechanisms in Sec. Ullhave to be compared 
in their relative strength for a given set of control pa- 
rameters. Here we focus on the phenomenological model 
from Hong et al. 24] and the approach based on kinetic 
theory from Trujillo et al. [23\ . 

We could not compare our results to predictions of 
hard sphere lattice models [5l|, [s^, [H, [13], which are 
based on a statistical mechanics approach to granular 
media [ssj . Their main control parameter is a configura- 
tional temperature, which we were unable to determine 
in our experiments and simulations. 



A. The model of Hong, Quinn and Luding 

While Hong et al. do not deny the presence of other 
segregation mechanisms in specific situations, they pro- 
pose that a theory based on void filhng (percolation in 
their nomenclature, discussed in Sec. Ill A 1[) and the con- 
densation mechanism (Sec. Ill C ip will provide the gen- 
eral frame work to predict if the BNE or the RBNE is 
observed. Interpreting the results from Rosato et al. [s^ . 
they assume void filling to be controlled by evf 



tnp camcta 



(4) 



Equating Cvf and Ccond (Eq- El) results in an equation 
for the phase boundary between the two effects which is 



PS 
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(5) 



with the RBNE if the left hand side is larger than the 
right hand side and the BNE for the opposite case. Hong 
et al. study however the more special case where the 
layer height ratio is chosen to be equal to the diameter 
ratio hi^/hs — di^/dg . Consequently they obtain for the 
phase boundary 



dL 
ds 



PS 

Pl' 



(6) 



which is in agreement with their molecular dynamics sim- 
ulations (23 | and a variational approach minimizing a bi- 
nary mixture free energy functional 56]. 

An experimental test of Eq. [5] or [5] is complicated by 
the fact that the condensation mechanism can only be ob- 
served for an intermediate granular temperature, which 
is difficult to achieve in experiment. Consequently Breu 
et al. Q did a first study of Eq. [5] by manually prepar- 
ing binary mixtures in the layer sequence predicted to be 
unstable for this diameter and density ratio. Then they 
searched for a combination of the control parameters / 
and A where the configuration changed to that predicted. 
This search succeeded for 82% of the 178 different particle 
combinations tested. However, they also observed that 
for the same mixture both BNE and RBNE can occur, 
depending on the driving parameters. Other experimen- 
tal [12, [531 and numerical ^28^ studies could not verify 
the existence of the condensation mechanism. 



i 



pLimp ~ 





m 



FIG. 1: Left: Experimental Setup. The particles are shaken 
vertically in an evacuated square glass tube. Pictures are 
taken from the top and the bottom (using a mirror) of the 
granular sample. Right: An image of the bottom plate. All 
large particles detected by the image processing are marked 
with black circles. 



of them are geometric in nature: the static compressive 
force (discussed in Sec. IH A2|) and the dynamic tensile 
force, which they compare with the void filling mech- 
anism (see Sec. IH A ip . The third term stems from the 
non-cquipartition of energy and is proportional to Tl/Ts; 
it is called pseudo-thermal buoyancy (Sec. Ill C 2p . Fi- 
nally, there is the classical Archimedean buoyancy, which 
is not relevant for particles of same material densities. 
The model of Trujillo et al. does not account for convec- 
tion or spatial inhomogeneities of the kinetic temperature 
in the sample. 

For the two-dimensional case they find that for small 
dh/ds the static compressive force, which leads to RBNE, 
is stronger than the other two terms, which favor BNE. 
However, the pseudo-thermal buoyancy grows faster with 
increasing dj^/ds, so there exists a critical diameter ratio 
above which the system exhibits BNE. They also predict 
that increasing the fraction of large particles will strongly 
increase the range of diameter ratios where the RBNE is 
observed. The same effect should apply for lowering the 
packing fractions by means of stronger shaking. 



IV. METHODS 



A. Experimental setup 



B. The model of Trujillo, Alam and Herrmann 

All segregation mechanisms discussed in this paper 
originate from the presence of a gravitational field paral- 
lel to the direction of shaking; in the framework of kinetic 
theory the influence of this field was first studied by Jenk- 
ins and Yoon [i^ . Starting from the same momentum 
balance equations, Trujillo et al. psj decompose the net 
gravitational force into four physical mechanisms. Two 



The experimental Setup is shown in Fig. [T] The par- 
ticles are shaken in a 25 mm x 25 mm square glass tube 
(Pyrex® Borosilicate) with a height of 360 mm. The bot- 
tom plate is made out of sapphire, the top cover of acrylic 
glass. The cell is mounted on an electromechanical shaker 
(Vibration Test Systems vglOOc) and shaken sinusoidally 
with an amplitude A and frequency /. The dimensionless 
shaking acceleration is F = {2Trf)'^A/g, where g is the ac- 
celeration due to gravity. F is controlled by a feedback 
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loop to better than 0.4%. 

All particles were made of brass (alloy 260) with a 
density p of 8.4g/cm'^. The diameter of the large 
particles was 2.38 mm; the small particle diameters dg 
were 0.79 mm, 1.19 mm, and 1.59 mm, corresponding to 
diameter ratios of 3, 2 and 1.5. 

All experiments were performed with the same proto- 
col: after filling with particles, the container was sealed 
and evacuated to a pressure below 8 Pa. Then a sequence 
of measurements was performed, each corresponding to 
a specific set of driving parameters / and F. Each in- 
dividual measurement started with the generation of an 
approximately mixed state by shaking 30 s at / = 15 Hz 
and r — 4.5. Then the particles were shaken with the 
specific / and F until a steady state was reached. After 
the shaker was stopped the particles came immediately 
to rest and images were taken from the top and bottom 
surface of the sample with a resolution of 40 pixels per 
di^. Three more sets of images were taken, each after 
shaking for another 20 s. The average number of parti- 
cles counted in the four sets of images was considered one 
independent measurement. 

To identify the position and type of all particles we 
use two images taken with different illumination. A first 
image, taken with a single point light source, marks the 
center of each particle with a light spot. A second image, 
taken with diffuse illumination, is then used to distin- 
guish between small and large particles. The algorithm 
makes this decision automatically by comparing the cor- 
relation between every found particle and template im- 
ages of small and large particles. The right side of Fig. [T] 
shows a result of this image processing. The error rate of 
the algorithm, as determined by comparing the number 
of large particles to the results of counting them by eye, 
is smaller than 6%. 

To quantify the amount of convection we took addi- 
tional images of a sidewall of the sample during the shak- 
ing process. Using a Phantom high speed camera we re- 
solve the particle dynamics with a resolution of 24 pixels 
per c/l at 1077 frames per second. The high frame rate al- 
lowed identification of the motion of individual particles. 
A small depth of field helped to distinguish particles near 
the sidewalls from particles in the bulk. Figure [5] gives 
four examples. 

In order to get reproducible results the brass particles 
had to be treated in two ways. First, all new particles 
were aged by shaking them for 24 hours at 20 Hz and F — 
5 to exclude effects due to work hardening. Second, after 
about every 8 hours of shaking the particles were cleaned 
for 10 minutes in an ultrasonic bath with a concentrated 
solution of Alconox powdered cleaner. This removed the 
fine metal dust which had accumulated in the sample. 
The resultant reproducibility is illustrated in the top row 
of Fig. [3] 




FIG. 2: Side view of the shaken sample while it is most ex- 
panded. The labels give the shaking frequency and amplitude, 
corresponding to F = 20, 1.7, 2.5, and 7, respectively. The 
container is filled with 128 large and approx. 1040 small brass 
spheres (corresponding to equal volumes); their diameter ra- 
tio is 2. The elongated appearance of the particles in the 
horizontal direction is a consequence of the illumination. The 
height of the images is one-eighth of the overall tube length; 
this height is large enough so that the spheres do not collide 
with the top of the container. 



B. Molecular dynamics simulation 

A molecular dynamics code described in [5^ is used 
for the simulations. In the simulation the container size, 
number of particles, diameter ratio, mass ratio, and shak- 
ing parameters are all identical to those in the experi- 
ment. Collisions are assumed to be binary and instan- 
taneous to permit the use of an event driven algorithm 
to decrease calculation time. The velocities before and 
after an event are calculated using the collision model 
proposed by Walton, which does not take into account 
the details of an interaction but does account for friction 
and rotation [59| . 

For both particle-particle and particle-wall collisions, 
three constants each (a total of 6) must be defined in 
this model. First is the coefficient of normal restitution, 
e, which is the ratio of the relative normal velocity (i;„) 
after the collision to the value before the collision. This 
coefficient in both cases varies with w„ as described in 
[ssl : the coefficient of restitution is the maximum of e and 
1 - (1 - e){vn/^/gdf/'^. The chan ges in relative surface 
velocity are determined by the coefficient of friction /i 
and the minimum coefficient of tangential restitution /3. 
The ratio of the value after to before the collision is the 
maximum of —j3 and /i times the normal impulse. This 
allows a transition from sliding to rolling contact. 



C. Comparison of experiment and Molecular 
dynamics simulation 

The collision parameters were systematically varied 
until optimal agreement was obtained between simula- 
tion and experiment; the results are shown in Fig. [3] lower 
row. This led to a choice of coefficient of normal restitu- 
tion, e = 0.78, and friction coefficient, ^ = 0.3, for both 
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FIG. 3: (Color online) Top row: Test of the reproducibility 
in three experimental runs, each starting with washed par- 
ticles as described in Sec. IIV Al Open symbols correspond 
to the average number of large particles detected at the bot- 
tom; closed symbols to the number on the top. Each set of 
points represents experiments with identical shaking accelera- 
tion. Bottom row: Comparison of average of the experimental 
results in the top row (squares) and the molecular dynamics 
simulations (triangles). The sample contains 64 large and 
1540 small brass spheres with a diameter ratio of 2. 



particle-particle and particle-wall collisions. 

For particle-particle interactions, the minimum tan- 
gential restitution, (3p, was set to 0.5; the results were 
only weakly dependent on Pp. Critical to the match was 
f3w, the tangential restitution for particle- wall interac- 
tions, which was set to 0.9. This sensitivity to (3^ hints 
towards the importance of sidewall-driven convection in 
the system. For a detailed discussion of the significance 
of (3 in this collision model, see [60j . 

To mimic the experiment, F was reduced to 0.5 after 
1000 cycles for simulations at 15 Hz and 40,000 cycles for 
80 Hz. In the event drive simulation the particles did not 
fully come to stop as quickly as they did in the experi- 
ment. The simulation continued until the time between 
collisions with the plate became too small to be handled 
by the algorithm. For F = 2, / = 15 Hz and F = 8, 
/ = 80 Hz, this collapse happened quickly as the layer 
was dense and close to the plate. The results from these 
simulations differ from experiment more than for higher 
shaking amplitudes as some of the dynamics involved in 
the layer coming to rest are not captured. 



The number of particles on top and bottom were de- 
termined in simulation from the heights of the particles. 
If the center a large particle was within O.Tc^l of the base, 
it was considered to be in the reverse state; if more than 
0.3(iL above the mean height of the small particles, in the 
Brazil nut state. This measure was determined by mak- 
ing images of the simulation results with the ray-tracing 
program Povray (sij . which were then analyzed in the 
same fashion as in the experiment. The cutoff values 
used minimized the difference between the two methods. 



V. RESULTS I: INFLUENCE OF DRIVING 
FREQUENCY AND AMPLITUDE 

In this section we will discuss the influence of the driv- 
ing parameters on the segregation in a "typical" sample 
consisting of each one monolayer (measured in units of 
c^l) of small and large particles with a diameter ratio of 
two: di^/ds = 2, Vl/Vs = 1, and /itotai = 2 c^l. We will 
first describe our experimental results and then use our 
molecular dynamics simulations to understand the role of 
convection and granular temperature inside the sample. 
In Sec. IV El we will then explain our findings in terms of 
the segregation mechanism discussed in Sec. [Hi Section 
IVII examines how these results change with the diameter 
ratio, the fraction of large particles and the overall filling 
height. 



A. Relaxation into a steady state 

An assumption of the measurement protocol described 
in Sec. IIV Al ls that the initial shaking period of 1 min is 
long enough for the system to reach a steady state. This 
is demonstrated in Fig. 2] where the shaker was stopped 
every 30 s to take images and count the number of large 
particles visible at the bottom and top surface. The time 
the sample needs to relax to a steady state decreases 
strongly with increasing acceleration. In experiments 
with F > 2.5 at 15 Hz and F > 8 at 80 Hz we reach 
a steady state in less than 1 minute, for lower accelera- 
tions we increased the shaking time accordingly. We also 
find this steady state in agreement with [J 0] to be his- 
tory independent: Starting from samples prepared in a 
BNE or RBNE configuration results in the same average 
number of particles. 



B. Dependence on the driving amplitude 

Figure [5] demonstrates that the segregation behavior is 
relatively complicated when plotted as a function of the 
shaking acceleration. For / = 15 Hz there is a crossover 
from a strong BNE to a weak RBNE at F w 4.5. In 
contrast for higher values of / we find the steady state 
always to be BNE and largely independent of F. 
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FIG. 4: Relaxation into a steady state exhibiting the BNE. 
The shaking frequency in the experiment is 15 Hz. The con- 
tainer is filled with 128 large and approximately 1040 small 
brass spheres (corresponding to equal volumes); their diame- 
ter ratio is 2. 



To characterize the segregation behavior we introduce 
the probability that an observed large particle is at the 
top surface, p = n-t / (nt+rii^) , where rit and rib correspond 
to the average number of large particles detected at the 
top and bottom surfaces. Consequently p = 1 identifies 
a complete BNE, p = a, RBNE state. 

The results for the probability p as a function of the 
driving amplitude A are shown in Fig.[6l with increasing 
shaking amplitude the system passes from a strong BNE 
to a weak RBNE. Qualitatively similar results have been 
found in the experiments of Breu et al. Q and the soft- 
core molecular dynamics simulations of Pica Ciamarra et 

al. m. 

Plots of p as a function of F (compare Fig. [5]), the 
average plate velocity or the particle velocity at lift off 
do not collapse the data as well as the dependence on A 
(Fig. [5]) ; therefore for the remainder of this paper we will 
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FIG. 5: (Color online) Experimental segregation results as a 
function of the shaking acceleration. The different symbols 
correspond to different shaking frequencies; the open version 
of a symbol represents the average number of large particles 
visible at the bottom, the closed version the number at the 
top. These experiments exhibit the BNE except for those 
with / = 15 Hz and F > 4. The container is filled with 128 
large and approx. 1040 small brass spheres (corresponding to 
equal volumes); their diameter ratio is 2. All data points 
for r > 2.5 are averaged over 19 independent experiments. 
For the measurements at 1.7 < F < 2.5 we used relaxation 
curves of the type shown in Fig. |4l where the initial shaking 
time necessary to reach a steady state (up to 15 minutes) was 
discarded. 
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FIG. 6: (Color online) Experimental segregation results as a 
function of the shaking amplitude. The data presented are 
the same as in Fig. [5l (pmix = 0.523 ^63]) 



describe our results using A. 



C. The role of convection 

Convection driven by the sidewalk of the container 
cannot be switched off in experiment. To estimate 
the relative importance of this mechanism we therefore 
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FIG. 7: (Color online) At 15 Hz (to the right of the vertical 
dotted line): simulations with frictional sidewalls (o) agree 
with experiment (•), while simulations with periodic bound- 
ary conditions (□) differ. This shows the importance of side 
wall driven convection rolls at this frequency. At 80 Hz (to the 
left of the vertical dotted line): the simulations are indepen- 
dent of the boundary conditions. The experimental results 
are the same as in Fig. [S] All simulation results are averaged 
over 10 runs starting from different initial conditions. 



rely on the molecular dynamics simulations described in 
Sec. IIVBI where we can suppress convection by using 
periodic boundary conditions for the lateral walls of the 
container 28]. Figure [7] shows again good agreement be- 
tween experiment and simulations with frictional side- 
walls for all but the smallest shaking amplitude at both 
15 Hz and 80 Hz. However, the role of convection is 
different at the two frequencies. 

Switching off convection in the simulations at 15 Hz 
results in a significantly lower probability for a large 
particle to be on top. For shaking amplitudes smaller 
than 1.2 dj^, switching off convection makes the system 
go from BNE to a mixed state. At higher amplitudes 
another RBNE mechanism eventually becomes stronger 
than convection so that the overall effect even with con- 
vection becomes RBNE. 

Simulation results at 80 Hz with frictional sidewalls 
and periodic boundaries agree quite well, showing that 
convection is not responsible for the strong BNE observed 
there. This agrees with the experimental observations: 
in high-speed movies of the top and bottom surface of 
the sample the driving mechanism was identified as void 
filling. 

This dependence of the segregation mechanism on 
/ agrees with the previous experimental observations: 
Vanel et al. [l^l found the motion of an intruder to 
be governed by convection for frequencies below 40 Hz 
and by void filling above. Other experiments observ- 
ing convection were performed with well separated taps 
[la nil, HO, EH, [13 or at 10 Hz 

Figure [S] shows the experimentally measured average 
vertical velocities in the layer adjacent to the sidewalls. 




1 2 3 

shaking amplitude in large particle diameters A/d|_ 

FIG. 8: Convection becomes important as the shaking ampli- 
tude increases. These measurements, performed using high- 
speed images as in Fig. [51 show the downward velocity of the 
particles at the sidewalls of the container. The sample is the 
same as in Fig. [S] The velocities are averaged over all phases 
of 6 cycles (15 Hz) or 30 cycles (80 Hz). 



For small A the observed velocities are comparable in size 
to the measurement errors. However, for / = 15 Hz and 
A > di,, the small particles start to become dragged down 
significantly, which is in agreement with the findings in 
Fig. El For shaking amplitudes larger than 1.5 di^, large 
particles also get entrained in the downward fiow. At 
these amplitudes the system shows a slight RBNE, so we 
are unable to decide if convection still drives the system 
towards BNE or if it has become a mixing mechanism. 



D. The granular temperatures of the two species 

Figure [9] shows the height dependence of the number 
density and the granular temperature in simulations with 
periodic boundaries. Because gravity and shaking differ- 
entiates between vertical and horizontal motion, we an- 
alyze the granular temperature in these two directions 
separately. The upper frame characterizes the fields at / 
— 80 Hz and A — 0.33 d^; their behavior is independent 
of the phase of the driving. The center of mass of the 
large particles is higher than that of the small particles, 
corresponding to the BNE found after switching off the 
shaking. While this BNE is generally attributed to the 
void-filling mechanism, it is interesting that the center 
of mass of the large particle has also the same height as 
the minimum of granular temperature in both horizontal 
and vertical directions. This hints that the realm of ther- 
mal diffusion and therefore kinetic theory might extend 
to relatively dense samples. 

At / = 15 Hz and A ~ 3.25 the results do depend on 
the phase of the shaking cycle. However, after the initial 
shock wave induced by the plate hitting the sample has 
traveled upwards fs^l, the remaining 5/6 of the shaking 
cycle is well represented by the data in Fig. [HI Two fea- 
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FIG. 9: (Color online) Number density and granular tem- 
perature in vertical and horizontal directions in simulations 
with periodic boundaries. As predicted by thermal difltusion, 
the larger particles accumulate preferentially at the minimum 
of the granular temperature. The sample contains one mono 
layer (measured in units of di,) each of small (o) and large (a) 
particles with a diameter ratio of two. The RBNE observed 
in the experiment at 15 Hz is in accord with the relative posi- 
tions of the arrows, which show the height of the the centers 
of mass for the small and large particles. The granular tem- 
peratures were only computed for heights where there was in 
average at least one particle in a height intervall of of the 
given type at a time; their values are scaled by the energy 
needed to raise a large particle its own diameter. Zero height 
corresponds to the lowest position of the vibrating plate. Data 
are averaged at a single phase of the cycle for 1200 (80 Hz) 
or 434 (15 Hz) cycles. At 15 Hz this phase corresponds to the 
most expanded sample at the upper turning point; at 80 Hz 
no dependence on the phase was observed. 



tures are characteristic for this "free flight" phase: first 
the center of mass of the large particles is always below 
that of the small particles; this agrees with the RBNE 
found in experiment and simulation. Second, the cen- 
ter of mass of the large particles is at approximately the 
same height as the minimum of granular temperature in 
both horizontal and vertical directions. 

These two features point clearly to the thermal diffu- 
sion mechanism suggested by Hsiau and Hunt ^45,] to be 
the cause for the observed RBNE. Minima of granular 
temperature arc a generic feature of vertically vibrated 
granular samples 37, 38, 39, 43], so thermal diffusion can 
be expected to be a relevant segregation mechanism in 
many systems. 

Characterizing the degree of non-equipartition is dif- 
ficult due to the strong spatio-temporal dependencies of 
the granular temperatures of both species. For a first 



estimate we have computed a weighted ratio: 

Ts ^ Ts{z) 

where ni^{z) is the height dependent number density of 
large particles. Averaged over the whole cycle, the sim- 
ulations at / = 15 Hz and A = 3.25 c^l yield a temper- 
ature ratio of 5.4 ± 2.2 in vertical direction and 2.8 ± 
0.6 in horizontal direction. The larger standard devia- 
tion in vertical direction results from a threefold increase 
of Tl/Ts during the initial shock wave. At / = 80 Hz 
and A = 0.33 cZl the temperature ratio is 2.2 in vertical 
direction and 1.7 in horizontal direction. 

Predictions of Tl/Ts from kinetic theory assume a di- 
lute granular gas either driven by a stochastic thermostat 
(Barrat and Trizac ^5^) or shaken at high frequencies and 
high amplitudes (Brcy et al. [48]). Figure [2] demonstrates 
that a granular gas is most closely realized by shaking 
with / = 15 Hz and A = 3.25 di^. With our diameter 
ratio, coefficient of restitution, and number of particles, 
the polynomial equation given by Barrat and Trizac pre- 
dicts Tl/Ts = 2.62, which is in good agreement with our 
result of 2.8 ± 0.6 for the horizontal temperature ratio 
(which depends less on the phase of the shaking cycle). 
The cubic equation for the ratio of the mean square ve- 
locities given by Brey et al. [11] results in Tl/Ts = 1.28; 
however, their result was derived for the case of a sin- 
gle large intruder, not equal volumes of small and large 
particles. We will return to this question in Sec. IVI Al 

E. Relevance of the different mechanisms 

The following overall picture emerges: at high frequen- 
cies and small shaking amplitudes the sample shows a 
strong BNE due to void filling. At low frequencies and 
small shaking amplitudes, sidcwall driven convection 
also leads to a BNE. However, with increasing shaking 
amplitude, thermal diff'usion becomes relevant: be- 
cause the minimum of granular temperature is lower than 
the center of mass of the whole sample, thermal diffusion 
results in a RBNE that eventually becomes stronger than 
the BNE due to convection. 

By Eqs.OandlH the control parameters for both con- 
densation and void filling have the same numerical value 
of 8. Therefore, the theory of Hong et al. predicts a 
strong BNE except for some intermediate driving range 
where condensation takes place and is of equal strength 
as void filling, which should lead to a mixed state. Fig- 
ure shows that a condensation of large particles is only 
conceivable for shaking amplitudes smaller than 1.16 cZl- 
Neither the experimental data in Fig. [5] nor the simula- 
tions in Fig. [7] show signs of condensation taking place at 
any frequency and amplitude. 

The effect of non-equipartition and static com- 
pressive force on the RBNE at high shaking amplitudes 
can be more appropriately studied by their dependence 
on the diameter ratio, as described in the next section. 
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VI. RESULTS II: INFLUENCE OF DIAMETER 
RATIO, VOLUME RATIO, AND TOTAL LAYER 
HEIGHT 

A. Dependence on diameter ratio 

The dependence of the experimental segregation re- 
suhs on the diameter ratio is shown in Fig. [TOl At / = 
80 Hz a change of the diameter ratio from 1.5 to 3 al- 
ways increases the BNE, which is in agreement with our 
interpretation of void filling as the relevant mechanism. 

At / = 15 Hz the situation is more complex: for A = 
1.16 c^L, the system exhibits BNE; its strength increases 
with diameter ratio. At larger shaking amplitudes we 
observe RBNE at small and BNE at larger diameter ra- 
tios. This general behavior is the same for a single large 
intruder and equal volumes of large and small particles; 
however, the effects are more pronounced for the single 
intruder. 

We know from Fig. [7] and [H that the results at / = 15 
Hz involve at least two mechanisms: convection and ther- 
mal diffusion. Since we want to use the dependence on 
the diameter ratio to study the role of non-equipartition 
and static compressive force, we first try to isolate the 
contribution of convection. In analogy to Fig. [7l we do 
this by comparing the experimental results with simu- 
lations with rigid and periodic boundaries, as shown in 
Fig. [TT] We consider only the case of the largest shaking 
amplitude because it can be expected to resemble best 
the driving assumed in [2^, . 

The agreement in Fig. [TT] between the experimental re- 
sults and the simulations with rigid boundaries decreases 
for diameter ratios of 1.5 and 3. This is probably due to 
two reasons: First we kept the collision parameters in 
the simulation at the values obtained from the fit with 
data for a diameter ratio of two; however, the coefficient 
of restitution is known to change with the particle di- 
ameter [ssf. Second, the large size discrepancy leads to 
increased probability of multi-particle collisions, thus the 
Tc model of Luding and McNamara was implemented for 
the simulations with a diameter ratio of 3 In this 

scheme, the restitution coefficient was set to 1 if a particle 
encountered a second collision within a short time frame. 
Still only a subset of the simulations ran successfully. 

Despite these caveats Fig. [TT] shows that the BNE for 



TABLE II: Comparison of the ratio of granular temperatures 
Tl/Ts between our simulations and the kinetic theory based 
predictions by Barrat and Trizac [50l |. The samples consist of 
1 layer of each of small and large particles and are shaken at 
/ = 15 Hz and A = 3.25 d^. The simulations use periodic 
boundary conditions; Tl/Ts is calculated using Eq.[7]and the 
horizontal temperatures. 





1.5 


2 


3 


Barrat and Trizac 


1.51 


2.62 


9.11 


MD simulation 


1.3 ± 0.1 


2.9 ± 0.6 


8.9 ± 3.8 



a diameter ratio of 3 is due to an increase in convection. 
This effect can also be seen in the average particle veloc- 
ity measured at the sidewall of the container in Fig. [12] 
At small shaking amplitudes the convection is again neg- 
ligible as in Fig. [8] For A > di^ the difference in the 
downward velocity for small and large particles increases 
significantly with the diameter ratio. 

Aside from the effect of convection, the main result 
in Fig. [TT] is that the slight RBNE in strongly shaken 
mixtures is relatively independent of the diameter ratio. 
Yet non-equipartition increases by a factor of six for the 
diameter ratios studied, as shown in Table [IT] which also 
shows that Tj^/Tg is in good agreement with the kinetic 
theory based predictions of Barrat and Trizac [l^l . 

Together this relative independence of p on c?l /dg and 
the strong dependence of Tj^/Tc^ on di^/dg indicates that 
non-equipartition has only a weak influence on our re- 
sults. This is in contrast to the predictions of Trujillo et 
al. [25j|, who expect the static compressive force which 
leads to RBNE to be dominant at small diameter ratios. 
However, non-equipartition grows faster with increasing 
diameter ratio than the static compressive force, so the 
system should show a crossover to BNE. 

For the reasons given above, we were unable to gather 
meaningful statistics for simulations of a single intruder 
in a container with periodic boundaries. Therefore, we 
cannot comment on the predictions of Brey et al. [48|, 
who found for an intruder that non-equipartition leads 
to a RBNE that increases in strength with increasing 
diameter ratio. 



B. Dependence on the volume ratio Vl/Vs 

The dependence of the segregation results on the com- 
position of the sample measured as the ratio between the 
volumes of the large and small particle Vl/Vs is shown 
in Fig. 1131 For each volume ratio the same behavior as 
in Fig. [6] can be observed: the probability p of a large 
particle to be at the top of the sample decreases with in- 
creasing A. For volume ratios in the range from a single 
intruder up to 20% large particles, p depends only weakly 
on Vl/Vs. 

At higher volume ratios, p shows a trend from BNE to- 
wards a mixed state with increasing Vl/Vs. This effect is 
in agreement with the suggested segregation mechanism: 
at / = 15 Hz the size of the convective downstream layer 
increases and more large particles can get entrained [3l| , 
while at / = 80 Hz an increase of Vl/Vs results in larger 
voids in the sample and therefore the void-filling mech- 
anism becomes less effective. A comparison of the two 
frequencies shows that convection is more impaired by 
an increase in Vl/Vs than void filling. 

An increasing RBNE with increasing diameter ratio 
was not found for any of the combinations of volume ra- 
tios and driving parameters in Fig. 1131 Therefore, a static 
compressive force is not needed to explain our results. 

Figure [T^] provides another test of the theory of Hong 
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FIG. 10: (Color online) Dependence of the observed segregation on the diameter ratio. The samples in the upper plot contain 
a single large particle in 2 layers (measured in di,) of small particles. The samples in the lower plot consist of equal volumes of 
each 1 monolayer of small and large particles. All data points are averaged over 19 independent experiments. 



et al. [131 . According to Eq. [5] the condensation mecha- 
nism leading to RBNE will become stronger than the void 
filling mechanism for volume ratios larger than 1. The 
right-most points in Fig. [T3]have Vl/Vs = 1.5; therefore, 
we should see here at some intermediate driving regime 
a clear drop of from a BNE to RBNE and back to BNE 
for a further increase in A. 

While no condensation was found at an intermediate 
driving range, this leaves open the possibility of conden- 
sation at even smaller driving amplitudes. Hong et al. 
consider a particle as being fluidized if it has the ability 
to exchange its position with its neighbors. Therefore, 
the mean square displacement of particles in a condensed 
phase should show a plateau at a small distance. Figure 
1141 shows the development of the mean square displace- 
ment of large and small particles in simulations with pe- 
riodic boundary conditions. Only at a frequency of 80 
Hz and an amplitude of 0.07 c^l (which pushes the lim- 
its of our simulation technique as discussed in Sec. IIVC|) 
do the large particles seem to be caged by their neigh- 
bors, while the small particles are still mobile. Accord- 
ing to the condensation mechanism, this should result in 
a RBNE. However, in experiment we find a clear BNE 
with p = 0.67 (at a Vi,/Vs of 1.5). This steady state is 



reached after approximately 20 minutes (corresponding 
to 10^ shaking cycles) and therefore unlikely to be found 
in simulations. Similar experiments at / = 15 Hz and 
shaking amplitudes down to 0.61 c^l do not yield any 
signature of the condensation mechanism. 



C. Dependence on the total layer height 

The dependence of the segregation on the layer depth 
/itotai is shown in Fig. [151 For a diameter ratio of 1.5 
the number of large particles at the top is mostly inde- 
pendent of /itotai while the number at the bottom drops 
significantly with increasing ft-totai- At a diameter ratio 
of 3 there is a slight increase of the number at the top 
while the number at the bottom stays relatively constant 
at low levels. 

The absence of the RBNE in deeper layers was also 
noted by Breu et al. It is compatible with thermal 
diffusion as the responsible mechanism for the RBNE in 
shallow layers, assuming that the relative height of the 
temperature minimum stays constant when /itotai is in- 
creased. Then in deeper samples the accumulation of 
large particles would no longer be visible at the bottom. 
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FIG. 11: (Color online) Dependence of the segregation on 
the diameter ratio with and without convection. The closed 
circles (•) correspond to the experimental results for f — 15 
Hz, A — 3.25 dL and equal volumes of small and large parti- 
cles(Fig. llOfl . The open circles (o) are the results of the MD 
simulations using frictional sidewalls. The open squares (□) 
represent the simulation results if periodic boundary condi- 
tions are used to suppress convection. All simulation results 
are averaged over 10 runs starting from different initial con- 
ditions. 




12 3 

shaking amplitude in large particle diameters 

FIG. 12: (Color online) Segregation due to sidewall driven 
convection increases with diameter ratio. The mean vertical 
velocity for small and large particles was measured in exper- 
iment at the sidewall of the container for diameter ratios of 
1.5 (T) and 3 (•). Small and large symbols represent the 
respective particles. Both samples contain equal volumes of 
small and large particles. The velocities are averaged over all 
phases of 6 cycles (15 Hz) or 30 cycles (80 Hz.) 



The condensation mechanism does not explain the van- 
ishing of the RBNE in deeper layers as the condensed 
fraction of large particles is assumed to accumulate at 
the bottom independent of /itotai- At a diameter ratio of 
3, void filling and convection are much stronger contribu- 
tions to the final result; therefore, the relative infiuence 



of /itotai is smaller. 

Note that the transition from two to four layers in- 
creases the relaxation times. Though we increased the 
initial shaking time from 1 to 3 minutes, experiments 
with / = 80 Hz, A = 0.13 c^l still had not reached their 
steady state and were therefore not included. 



VII. CONCLUSIONS 

We have studied the segregation of vertically shaken 
binary mixtures of brass spheres in experiment and event- 
driven simulation. By using the same material for small 
and large particles and evacuating the sample container, 
we eliminated Archimedean buoyancy, friction, and pres- 
sure effects from the list of potential segregation mecha- 
nisms. We have found that three mechanisms are suffi- 
cient to explain our results: 

(1) At a shaking frequency of 80 Hz all samples ex- 
hibited a strong Brazil nut effect due to the geometri- 
cal mechanism called void filling. The strength of the 
segregation increases with diameter ratio and decreases 
weakly with the shaking amplitude. Increasing the num- 
ber of large particles from a single intruder up to a mix- 
ture of equal volumes of large and small particles results 
only in a slight decrease of the strength of the effect. 

(2) At a shaking frequency of 15 Hz and shaking ampli- 
tudes up to 1.5 large particle diameters, a strong Brazil 
nut effect is induced by side wall driven convection. 
The strength of the effect increases strongly with diame- 
ter ratio and decreases with the relative amount of large 
particles. 

(3) If the shaking amplitude is increased above 1.5 
large particle diameters, the sample starts to resemble 
a granular gas and the thermal diffusion mechanism 
becomes significant. Thermal diffusion describes the ten- 
dency of large particles to accumulate in the minimum 
of the granular temperature profile. In a shallow layer 
of height two large particle diameters and for diameter 
ratios smaller than three, thermal diffusion eventually 
becomes stronger than convection, which leads to the re- 
verse Brazil nut effect. Thermal diffusion seems to be 
only weakly dependent on the diameter ratio. If the to- 
tal layer height is increased to four large particle diame- 
ters, the temperature minimum is within the sample and 
no increased number of large particles at the bottom is 
observed. 

While we probed the parameter regimes where the 
static compressive force and the condensation 

mechanisms were expected to be relevant, we did not 
find any signature of these two segregation mechanisms. 
The apparent agreement between the condensation-based 
model in [2^ and the experimental phase diagram sep- 
arating BNE and RBNE in Q can be explained by the 
fact that they changed both frequency and amplitude of 
the driving deliberately. Then both BNE and RBNE can 
occur in the same sample, as shown above and also in 
We also find non-equipartition to have no discernible 
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FIG. 13: (Color online) Dependence of the observed segregation on the volume ratio between small and large particles. All 
samples have a total volume of 2 layers (measured in cJl)- The left-most data points in each frame correspond to a single large 
intruder. All results are averaged over at least 18 independent experiments. 



influence on our results. However, this might be due to 
either the relatively high coefficient of restitution (0.7) 
or to the fact that we are not shaking vigorously enough 
for those effects to be observed. 

The theories presented in [23, UHl lack two of the im- 
portant mechanisms (convection, thermal diffusion) and 
at the same time they include mechanisms which are 
not relevant (static compressive force, non-equipartition, 
condensation); hence these theories can not explain our 
experimental and numerical results. 

In the limit of strong shaking, we believe that kinetic 
theory-based approaches have the potential to explain 
segregation, though the inclusion of convection might 
prove daunting. A good starting point for a microscopic 



theory of void filling might be the "random fluctuating 
sieve" theory of inclined chute flow [6^, as pointed out 
in [31]. 

Note added in proof. The assumption in VI C that the 
relative height of the temperature minimum is indepen- 
dent of the total number of spheres is supported by the 
analysis in Ref 64]. 
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FIG. 14: (Color online) Driving parameters adequate for the 
condensation mechanism were searched in simulations (with 
periodic boundary conditions) that yielded this plot of the 
mean square displacement. The diameter ratio is 2; the sym- 
bol size corresponds to the particle type. The sample size is 
2 layers with a volume ratio of large to small spheres of 1.5. 
All data points are averaged over 200 time steps. 
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